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' (—i '[ Abstract 

46 years ago the quark model replaced the Sakata model as the standard ex- 
planation of the hadron structure. The major alleged defect of the Sakata model 
(3jT)' was its prediction of just too many types of particles, which have not been seen in 

■ experiments. However, this allegation was made without detailed consideration 

of the forces acting between sakatons. In this article we suggest a set of pairwise 
sakaton-sakaton and sakaton-antisakaton potentials that describe stability and 
masses of strongly interacting elementary particles in a good agreement with 
Q_|. observations. 

1 The Sakata model 

oo , 
m . 

Today it is universally accepted that hadrons are made of quarks (q = u, d, s, c, . . .). 
The quark model forms the basis of quantum chromodynamics (QCD), which aspires to 
explain the nature of strong interactions. Almost all compound particles predicted by 
O \ the quark model have been found in experiments. Moreover, all observed particles have 

natural quark assignments: Mesons are quark-antiquark bound state (qq), and baryons 
are bound states of three quarks (gc/qOlj In spite of their well-known achievements, 
^ . the quark model and QCD have some questionable features. These theories make 

assumptions (fractional charges of quarks, color, gluons, confinement potentials, etc.), 
which cannot be directly observed and thus destined to remain suspect. Then it seems 
justified to explore other approaches to the explanation of hadrons masses and stability. 
One interesting proposal is the Sakata model [2] , which was rather popular before the 
"quark era". The Sakata model assumes that proton (p), neutron (n), A and A+ 



1 Some suggested tetraquark and pentaquark assignments are not universally accepted. All experi- 
mental data about hadron properties used in this paper were taken from [T] . Experimental values are 
emphasized by the bold font in this paper. 
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are the true elementary particlecl also called sakatons (a) [3J. To emphasize their 
similarity with quarks, we will denote the four fundamental sakatons by capital letters 
U, D, S, CE Each sakaton has its corresponding antisakaton (U, D, S, C) with the same 
mass and spin and opposite values of the electric charge, baryon charge, strangeness, 
and charm. 



Table 1: Stable baryons with their quark and sakaton structures. 



baryon 


Quark 


sakaton 


Exp. mass 




structure 


structure 


MeV/c 2 


n 


udd 


D (down) 


938 


P 


uud 


U (up) 


940 


A 


sud 


S (strange) 


1116 


A c + 


cud 


C (charmed) 


2285 


£~ 


sdd 


SDU 


1197 


S° 


sud 


SNN 


1193 


£+ 


suu 


SUD 


1189 




ssd 


SSU 


1322 


wo 


ssu 


SSD 


1315 


wo 

1—1 c 


csd 


CSU 


2471 


1— 'c 


csu 


CSD 


2468 


' — 'cc 


ccd 


ecu 


3519 


' — 'cc 


ecu 


CCD 


not seen 


n- 


sss 


SSSUD 


1672 


n° 


CSS 


CSSUD 


2698 


n£ 


CCS 


CCSUD 


not seen 


"cCC 


ccc 


CCCUD 


not seen 



The Sakata model assumes that sakatons interact with each other via short-range 
(few femtometers) potentials. All non-elementary hadrons are bound states of two 
or more sakatons. Various possible combinations are summarized in Table [2J Nuclei 
are composed of U and D sakatons (protons and neutrons). Mesons are sakaton- 
antisakaton (era) bound states. Compound baryons are sakaton-sakaton-antisakaton 
(excrex) or pentasakaton (aaaaa) bound states. Here we are interested only in baryons, 
which are stable with respect to strong decays. All of them are listed in Table [TJ Their 
decays are caused by flavor-changing weak interactions and their masses are lower 
than the sums of masses of constituents. One example of an unstable baryon state 

2 In this paper we will not discuss bottom and top particles, because full experimental picture is 
still lacking in those sectors. 

3 See Table [TJ We use symbol iV do denote collectively U and D sakatons. For example, NNS 
means either UUS or DDS. 



2 



omitted in Tabled] is the A ++ (= UUD) particle whose mass is 1232 MeV/c 2 . This is 
higher than the sum of masses of dissociation products p(940) + 7r + (140)(= U + UD). 
Therefore A ++ is a metastable resonant state in our model. The calculation method 
adopted in this work (see section [2]) can deal only with true bound states, therefore we 
will not discuss the A ++ and other resonances. 

The Sakata model avoids some problems characteristic for the quark model. The 
fundamental constituents of the Sakata model - the sakatons - are readily observable 
as normal baryons with integer charges, so there is no need for additional assumptions 
about "confinement". There is also no need to introduce "hidden" degrees of freedom, 
such as color and gluons. The short-range character of sakaton potentials means that 
strong interactions satisfy the important property of cluster-separability [I] , similar to 
electromagnetic and gravitational forces. 



Table 2: Bound states of sakatons 


sakaton content 


particle type 


examples 


antiparticle 


a 


baryon 


p, n, A , A+ 


a 


NN' . . . N" 


nucleus 


deuteron(= UD) 


nn' ...n" 


aa 


meson 


K+{= US) 


aa 


aaa 


baryon 


£-(= SDU) 


aaa 


aaaa 


tetrasakaton 


unstable? 


aaaa 


aaaa 


tetrasakaton 


unstable? 


aaaa 


aaaaa 


baryon 


fi-(= SSSUD) 


aaaaa 



The biggest problem of the Sakata model is that it seemingly predicts more types 
of particles than actually observed. Certain aaa combinations, which look acceptable 
from the point of view of the Sakata model, have not been seen in experiments. This 
refers, for example to NNS and UDS baryons with strangeness Furthermore, 
the simplest sakaton assignment of the Q~ baryon (baryon number = 1, charge = -1, 
strangeness = -3) is in the form of a pentasakaton f2~(= SSSUD). Then, from the 
principle of isotopic invariance, it seems that analogs of the Q~ particle should also 
exist, such as Q (= SSSUU) and ft°(— SSSDD). Why haven't they been seen in 
experiments? 

In order to answer these and other questions, it is important to have a realistic 
model of interactions between sakatons. The goal of this paper is to suggest an ap- 
proximate set of pairwise sakaton-sakaton and sakaton-antisakaton potentials and to 
calculate masses of their bound states - mesons and baryons. 

4 Reports about discovery of the exotic baryon 6 + (= UDS) are not credible pQ. 
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2 Computational model and results 



Matumoto and co-authors established [5j El [7J EJ |9] that masses of hadrons can be 
roughly calculated from the assumption of strong attraction in sakaton-antisakatons 
pairs (i.e., one a — a bond contributes about 1275-1740 MeV to the binding energy) 
and equally strong sakaton-sakaton and antisakaton-antisakaton repulsions. Binding 
energies of mesons are very high (above lGeV), because only the a — a attraction 
contributes there. Much lower binding energies are expected in 3-sakaton aaa and in 
pentasakaton aaaaa baryons. In the former case two attractive interactions o — a 
are balanced by one repulsion a — a. In the latter case there are 6 attractions vs. 
4 repulsions. Tetrasakatons aaaa are not likely to be stable because the number 
of repulsive and attractive pairs is equal in this case. Some instructive studies of 
multiparticle systems with pairwise interactions can be found in [101 E] . They suggest 
that stability of multi-sakaton states may depend on a delicate balance of masses of 
the constituents and shapes of their interaction potentials. 

The approximate non-relativistic Hamiltonian describing an AA-sakaton system can 
be written as 



where m i} p i} = |r$ — r 3 -| are masses and momenta of the sakatons and their relative 
distances, respectively. Interactions between sakatons were modeled as superpositions 
of two Yukawa potentials 



where Zi — +1 for sakatons and Z{ — — 1 for antisakatons. 

All calculations were performed using the stochastic variational method of Varga 
and Suzuki [121 H21 HI]- The FBS computer program was obtained from the CPC 
Program Library (Queen's University of Belfast, N. Ireland) and slightly modified to 
fit our needs. This program solves the non-relativistic stationary Schrodinger equation 
and yields accurate energies and wave functions of the ground and few excited states for 
systems of several (typically, 2-6) quantum particles interacting via pairwise potentials. 
Only states with the lowest total spin (s = for mesons and s = 1/2 for baryons) and 
zero orbital momentum were considered here. In all calculations masses of sakatons 



were fixed as m(N) = 940 MeV/c 2 @ m(S) = 1116 MeV/c 2 , m(C) = 2285 MeV/c 2 . 



5 The equality of masses of the U and D sakatons and the assumption that their interactions with 
other sakatons are the same (see Table UJ) imply that all calculated masses are invariant with respect 
to replacements in which all U sakatons are changed to D and all D sakatons are simultaneously 
changed to U. 




(1) 




(2) 
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Internally in the code these masses were expressed in units of the proton mass 940 
MeV/c 2 . Distances were measured in femtometers and energies in MeV. In this system 
of units h 2 /m = 41.47. The basis set selection procedure used iteration numbers 
Mq = 10, Kq = 50. Other computational parameters depended on the number of 
sakatons in the system as shown in TableEJ They were adjusted for the optimal balance 
between accuracy, convergence, and speed. The exact meaning of these parameters was 
explained in \1'2\ . 



Table 3: Computational parameters for the FBS code. b m i n /b max are mini- 
mum/maximum values of nonlinear parameters in Gaussian basis functions. 



Number of 


Basis set 




i 


Umax 


sakatons 


size, K 


(fm 


) 


(fm) 


2 


50 


io- 




10 


3 


250 


io- 




10 


4 


300 


10" 




10 


5 


500 


io- 




100 


Table 4: Optimized parameters of the potentials (T2 


). A&ndB 


are measured 


in MeV-fm; 


a and (3 are 


in fm 








Interaction 


A a 




B 




N — N 


617.8 0.091 




92.14 


0.359 


U-D 


570.2 0.091 




25.7 


0.094 


N-S 


530.0 0.108 




14.0 


0.49 


S-S 


446.7 0.118 




42.1 


0.444 


N-C 


397.5 0.102 




14.0 


0.49 


S-C 


340.8 0.12 




46.1 


0.444 


C-C 


317.0 0.118 




24.1 


0.484 



Our major goal is to optimize parameters A,a,B,(3 of the potentials ([2]). The 
optimization was performed in two steps. In the first step we fitted parameters relevant 
to interactions of U, D, S sakatons. The training set contained 24 species shown in 
Table [51 They included 4 ground states of mesons, 4 stable baryons, and 16 states, 
which are supposed to be unstable. The goal was to reproduce experimental masses 
of the 8 stable species as close as possible and, at the same time, do not allow the 
binding energy of the 16 unstable species to become positive. In the second step we 
froze the U — D — S parameters obtained above and varied interactions C — N, C — S, 
and C — C using the training set in Table [6j This set included C-containing particles: 
3 mesons, 3 stable charmed baryons, and 17 unstable species. The final optimized 
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Figure 1: Schematic shape of the U — D, U — D and U — D interaction potentials (full 
lines). The broken line shows a proposed modification of the U — D potential that can 
explain the bonding in nuclei. 

values of parameters A, a, B, (3 are given in Table HI Plots of the optimized U — D 
and U — D (same as U — D) potentials are shown in Fig. [TJ Potentials for other pairs 
of sakatons have qualitatively similar shapes. These interactions demonstrate rather 
strong attraction of a — a pairs and repulsion of o — a and a — a pairs in a qualitative 
agreement with Matumoto's guesses. 

The resulting masses of hadrons are shown in the third column of tables [5] and [6j 
The binding energies (B.E.) are in the 4th column and the lowest-energy dissociation 
products are in the 5th column. 

Ideally, the binding energies of unstable tetrasakatons and baryons^l must be equal 
to zero. In practice this can be achieved only with very large and diffuse basis sets, 
which allow the wave functions of dissociation products to separate widely, so that their 
repulsion is minimized. For computational reasons our basis sets were limited. This 
explains why some residual repulsion (reflected in negative binding energies from to 
-5 MeV) remained for several dissociated unstable species. Extremely large negative 
binding energies of UDUD, SDUD, SDUUD, CSUUD and CSUUD are explained 
by the fact that they have converged to metastable dissociated configurations (=local 
minima) UU+DD, SU+DD, SDD + UU, CSD + UU, and C+SU+UD, respectively. 

The next step is to consider properties that have not been used directly in the fitting. 

6 They are shown in the lower portions of Tables [5] and [S] 
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Table 5: Training set for no-charm particles. 



Particle 


Sakaton 


Mass (MeV/c 2 ) 


B.E. (MeV) 


Products 




structure 


calc./exp. 


calc./exp. 


calc./exp. 


7T° 


DD 


238/135 


1642/1745 


n + n/n + n 


7T + 


UD 


142/141 


1738/1739 


p + n/p + n 


K~ 


SU 


364/494 


1692/1562 


p + A°/p + A 


i] 


ss 


1095/548 


1137/1684 


A° + A°/A° + A° 


£+ 


SUD 


1210/1189 


48/67 


A° + 7r+/A° + 7r+ 


s° 


suu 


1260/1193 


44/58 


P + K-/A + 7T° 


wo 


SSD 


1314/1315 


166/295 


A + K°/A° + K 


Q- 


SSSUD 


1670/1672 


8/136 


S" + K / S- + K° 


unstable 


uuu 


1179 


-1/0 


P + 7T°/p + 7T° 


unstable 


UDU 


1082 


0/0 


P + 7T _ /p + 7T _ 


unstable 


uus 


1307 


-3/0 


p + K + /p + K+ 


unstable 


UDS 


1307 


-3/0 


p + K°/p + K° 


unstable 


sus 


1482 


-2/0 


A° + K+/A° + K+ 


unstable 


sss 


2211 


0/0 


A + T]/A° + 7] 


unstable 


UDUD 


483 


-199/0 


7T+ + 7T~ /7T° + 7T° 


unstable 


SDUD 


606 


-101/0 


K° +7T-/K- + 7T° 


unstable 


SDSD 


728 


0/0 


K° + K°/K° + K° 


unstable 


ssuu 


733 


-5/0 


K- + K-/K- + K- 


unstable 


SSUD 


731 


-3/0 


K~ + K°/K~ + K° 


unstable 


ssss 


2192 


-2/0 


7] + 7]/r] + 7] 


unstable 


UDDUU 


1223 


0/0 


P + 7T _ + 7T _ /p + 7T _ + 7T _ 


unstable 


SDUUD 


1499 


-147/0 


S+ + 7T-/S + 7T° 


unstable 


SSDUD 


1457 


-1/0 


S° + 7T-/S- + 7T° 


unstable 


sssuu 


1678 


0/0 


S- + A-/S- + K- 



First, we looked at the two charmed baryons whose existence is predicted by the quark 
model and whose experimental confirmation is still lacking. These are the fi+ and 
particles. We found that Q+(= CCSUD) is stable with the mass of 4430 MeV/c 2 
and binding energy of 25 MeV with respect to the fi+ -> S+ + iC°(= CC77 + STJ) 
dissociation channel. The calculated mass of f2^(= CCCUD) is 6099 MeV/c 2 , which 
means that this particle dissociates spontaneously as — >■ H++Z} + (= CCU + CD). 

Next we verified that all 133 possible tetrasakaton (crcraa) and baryon (aaa and 
aaaWa) species not presented in Tables [5] and |6] are unstable in our approach, as 
expected. 

Other interesting pieces of information, which have not been involved in the fitting, 
are the meson excitation energies. Note that the strongly attractive aa potential (see 
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Table 6: Training set for particles with charmed sakatons 



Particle 


Sakaton 


Mass (MeV/c 2 ) 


B.E. (MeV) 


Products 




structure 


calc./exp. 


calc./exp. 


calc./exp. 


D° 


CU 


2001/1897 


1224/1328 


p + A+/p + A+ 


Df 


CS 


2587/1968 


814/1433 


A° + A+/A° + A+ 


Vc 


CC 


3339/2980 


1231/1590 


A+ + A c + /A+ + A c + 


wo 

1 — 'c 


csu 


2606/2471 


43/308 


A+ + K-/A+ + K 


1 — 'cc 


ecu 


4092/3519 


194/636 


A+ + D°/A+ + D° 




CSSUD 


2852/2698 


118/268 


~ + + K-/S+ + K" 


unstable 


UDC 


2944 


-3/0 


p + D'/p + ~D- 


unstable 


uuc 


2945 


-4/0 


V + S°/p + D° 


unstable 


sue 


3121 


-4/0 


A +Tf/A° + D° 


11T1 cf Q nip 


ssc 

kJ v_y 


3705 


-2/0 


ao i n-/A° 4- D~ 


n n cf q nip 


CAW 


2427 


0/0 


A+ + 7T+ /A + 4- 7T + 


nnst^HIp 

U.J.J.Q U O) UlVj 


CUU 


2523 


0/0 


A + + 7T°/A + + 7T° 


unstable 


CSS 


3380 


0/0 


A+ + v/Af + r? 

c 1 lie 1 / 


unstable 


CSC 


4456 


-1/0 


A + r/ c /A° + r? c 


unstable 


CCS 


4872 


0/0 


A+ + D+/A+ + D+ 


unstable 


ccc 


5625 


-1/0 


A+ + r) c /A+ + rj c 


unstable 


cuus 


2365 


0/0 


D° + K+/D+ + K° 


unstable 


cu~ss 


2953 


-2/0 


D+ + K+/D+ + K+ 


unstable 


csss 


3685 


-3/0 


DJ + tj/DJ + jj 


unstable 


CUDUD 


2568 


-1/0 


A+ + 7T+ + 7I--/A+ + 7T° + 7T° 


unstable 


CSUUD 


2850 


-102/0 


2° + 7T+/S+ + 7T° 


unstable 


CSUDD 


2794 


-46/0 


H+ + 7T+/H+ + 7T+ 


unstable 


CSSDD 


2970 


0/0 


S+ + JfO/S+ + K° 



Fig. [TJ can accommodate a few stationary states that can be regarded as excitations of 
the ground-state meson. In Table [7] we show calculated masses of 3 lowest spherically 
symmetric (J p = 0~) meson states and compare them with experimental numbers 
where available. The same basis set was used for the ground and excited states. 
Obtained excitation energies of the order of several hundreds of MeV are roughly 
consistent with observed data. This gives us some confidence regarding the overall 
shape of the selected interaction potentials. 
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Table 7: Low mass states of some mesons with angular momentum quantum numbers 
S = L = J = 0. Masses are in MeV/c 2 . 



Sakaton 


1% 


mass 


2% 


mass 


3% 


mass 


structure 




calc./exp. 




calc./exp. 




calc./exp. 


DU 




142/140 


tt(1300) 


1480/1300 


tt(1800) 


1726/1816 


SU 


K~ 


364/494 




1669 




1909 


ss 


V 


1095/548 


77(1475) 


1974/1476 




2136 


cu 


D° 


2001/1897 




2942 




3117 


cs 


Df 


2587/1968 




3214 




3331 


cc 


Vc(lS) 


3339/2980 


Vc(2S) 


4282/3637 




4456 



3 Discussion 

The most important lesson of the above calculations is that the Sakata model is quali- 
tatively correct, at least in the part concerning masses of strongly interacting particles. 
With properly adjusted interaction potentials, this model correctly predicts the sta- 
bility of those species, which are found stable in nature. On the other hand, the 
unbound combinations of sakatons are exactly those, which were not seen in experi- 
ments. The calculated masses of stable particles (see Tables and E]) sometimes differ 
from experimental values by hundreds of MeV/c 2 . For example, masses of baryons are 
systematically overestimated. However, such discrepancies are expected due to our use 
of simplified 2-particle potentials (T5]). One can expect that true sakaton interactions 
have a more sophisticated form. 

For example, in our approach, U — D, U — U and D — D potentials are purely 
repulsive. This does not allow us to describe bound states like UD (deuteron) or 
UUDD (a-particle). It seems plausible that these interactions (especially the U — D 
potential) can be slightly modified so as to make them attractive at distances ~l-2 
fm (see broken line in Fig. [T]). Then it might be possible to reproduce the bonding 
of protons and neutrons in nuclei. Such a possibility is especially exciting as it would 
allow us to describe the stabilities of mesons, baryons and nuclei within the same set 
of sakaton interactions. 

Another missing piece is the absence of relativistic corrections that may include 
momentum-dependent, spin-orbit, spin-spin, and contact interactions. It is well-established 
that they can contribute up to several hundreds of MeV to the overall energy balance 
of mesons and nuclei. 

One can also add to (T5J) terms which change the number and/or types of particles. 
For example, terms like3 

7 Here u,u,d,d are annihilation operators for the U,U,D,D sakatons, and , , eft , d} are their 
creation operators. 
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V m i X oc v)v)dd + d)d uu + . . . (3) 

are responsible for the mixing of UU and DD states and for the mass splitting between 
7T° = l/y/2(UU-DD) and rj' « 1/V2(UU + DD) mesons PQ. Without interaction © 
particles 7r° and r/ have the same mass, while experimentally their masses are quite 
different: m(n ) = 135 MeV/c 2 , m(r/) =958 MeV/c 2 . This indicates the significant 
role of terms like fl3]). Generally, one can also expect the presence of interactions that 
lead to the mixings NN -h- SS -H- CC '. Our neglect of these interactions may partially 
explain the overestimation of masses of 7r°, 77, and rj c mesons. 

One may argue that Sakata's assumption of a fundamental point-like proton must 
be wrong because, being probed by truly point-like electrons, the proton demonstrates 
a sizeable charge radius of 0.877 fm. However, this experimental fact can be accommo- 
dated within the Sakata model as well. To achieve that, one can assume the presence 
of particle-number-changing interaction terms like 



Vunphys oc u ( u ( u [ u + d^d^u^u + u^uuu + v) ddu + . . . (4) 

in the Hamiltonian. In the classification of [15] these terms are called "unphys". If 
they are present, then single "bare" proton states u*\0) are not eigenstates of the total 
Hamiltonian. To make the theory sensible, one needs to perform a renormalization. If 
coefficient functions in the interaction © are properly chosenjl then all loop integrals 
are finite, the renormalization effects are finite too, and the "bare" proton becomes 
"dressed" by a cloud of virtual pairs and pions, thus acquiring a non-zero size [T6] . 

In spite of the deficiencies listed above, our results indicate a remarkable consistency 
between the quark and Sakata models: both models predict the same set of stable 
hadron states, § This suggests that Sakata's idea about the hadron structure has a 
non-vanishing fighting chance against the quark model. Further studies with more 
elaborate potentials would be certainly welcome. In addition to the masses of stable 
species considered here, these future studies should address resonances and scattering 
properties as well. 

I would like to thank Dr. Robert Wagner for reading the manuscript and for helpful 
critical comments. 
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